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EFFECTS  OF  RANDOM  SHADINGS,  PHASING  ERRORS, 
AND  ELEMENT  FAILURES  ON  THE  BEAM  PATTERNS 
OF  LINEAR  AND  PLANAR  ARRAYS 


INTRODUCTION 

Although  an  array  beamformer  with  known  element  positions  may  be  designed  for  good  sidelobes 
or  skirt  selectivity  behavior  by  choice  of  the  element  shading  (weights),  the  actual  array  response  will 
undergo  degradations  due  to,  e.g.,  element  position  movement,  delay  approximations,  element  gain 
quantization,  random  element  gains,  element  failure,  etc.  Here  we  will  investigate  the  effect  of  all 
these  random  perturbations  on  the  response  of  the  array  beamformer  for  a  single-frequency  plane- 
wave  arrival  and  for  individual  array  elements  with  omnidirectional  response.  Since  the  power 
response  to  a  single-frequency  plane  wave  is  itself  a  random  variable  at  each  angle  of  look,  we  will 
evaluate  its  mean  and  variance  as  a  function  of  the  look  angle;  the  element  locations;  the  statistics  of 
the  shading,  phasing,  and  failure  perturbations;  and  the  plane-wave  arrival  frequency,  propagation 
speed,  and  arrival  angle.  From  these  results  can  be  deduced  quantitative  tolerance  limits  on  the 
perturbations  in  order  to  realize  specified  sidelobe  levels. 

Some  previous  results  on  the  array  power  response  for  random  perturbations  of  the  element  gains 
alone  are  given  in  reference  1.  (This  reference  is  also  useful  for  additional  background,  motivation, 
and  interpretations.)  Then,  in  reference  2,  the  moments,  through  order  four,  of  a  sum  of  independent 
complex  random  variables  were  derived,  as  were  the  cumulants  through  order  six.  However,  both  of 
these  results  were  given  in  terms  of  the  moments  of  the  zero-mean  random  variables  of  each  com¬ 
ponent  in  the  sum;  this  form  is  rather  inconvenient  and  error  prone  when  calculating  array  per¬ 
formance.  Also,  the  beamformer  application  in  reference  2  was  limited  to  phasing  errors  only. 

Here  we  will  derive  the  moments  of  a  sum  of  independent  nonidentically  distributed  complex 
random  variables,  up  through  the  fourth  order,  in  complete  generality,  with  no  Gaussian  assump¬ 
tions.  Then  we  will  apply  these  results  to  both  linear  and  planar  arrays  and  give  examples  of  the 
performance  degradation  caused  by  perturbations  in  gain,  phasing,  and  element  failures.  Some 
results,  without  the  derivations  and  programs  contained  herein,  have  already  been  presented  in 
reference  3.  Additional  related  results,  which,  however,  do  not  cover  the  higher  order  moments 
considered  here,  are  presented  in  references  4  and  5. 

DEFINITION  OF  TERMS  AND  NOTATION 

Let  C  be  a  complex  random  variable.  The  average  value  of  C  is  denoted  here  by  two  equivalent 
notations, 

Me]  -  c,  (i) 

and  is  a  complex  quantity.  In  a  similar  manner,  we  have 

Av  { =  C3  ,  Mler]  «  Tcp\  (2) 


l 
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The  variance  of  complex  random  variable  C  is  defined  as 


V*rfc)  ■  |C-C|’  *  |cp-  |f|\ 


and  is  real.  Two  other  averages  of  interest  are 


Avf  (c-  O’]  .  (c-zj  -  ?- 
A» l icf}  *  W. 


The  final  variance  of  interest  is 


Vorfld']  *  =  |C\*  -  ]cf  . 


The  expression  in  (5)  requires  the  fourth-order  moment  in  (4)  and  constitutes  the  major  analytical 
problem  addressed  herein.  We  will  derive  expressions  for  all  the  quantities  in  (l)-(5). 

MOMENTS  OF  SUM  OF  RANDOM  VARIABLES 

The  particular  problem  of  interest  in  this  section  is  as  follows:  {z^f1  are  statistically  independent 
complex  random  variables,  which  are  not  necessarily  identically  distributed,  nor  are  they  assumed 
Gaussian.  A  sum  variable  C  is  defined  as 

c  *  <« 

and  is  complex.  We  wish  to  evaluate  the  various  averages  defined  previously  in  (l)-(5),  in  terms  of  the 
appropriate  moments  of  random  variables  {zk}.  These  moments  of  {zk}  are  presumed  known,  but 
they  need  take  no  special  form;  random  variable  zk  need  not  have  zero  mean,  for  example. 


We  have  immediately,  from  (6), 


Avfc]  *  C  - 


in  terms  of  the  means  of  { zk) .  Also  there  follows 

Av{c')  ■  c7  -  jf-KS  ■  fV  * 

•  +  (8) 

where  we  used  the  statistical  independence  of  zk  and  zn  for  k  #  n.  Continuing  on,  in  an  obvious 
fashion, 

MW] * ¥? ■  cc5 - 

■  i- (sy-  m1)  y  |  fxl*  *  f  (W-lsO  v  Icf, 
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where  we  used  (7).  From  (3)  and  (9),  there  follows 

Var|c]  *  ^  (  |%r  -  \%f)  =  ;  (10) 

and  from  (4),  (8),  and  (7), 

Av  ?  (C  -  ^T)3]  =  (11) 

An  alternative  expression  for  (9)  is  available  by  employing  (10): 

Av  l  |C|*]  *  |Av<C]|**  Vorfc}  (12) 

The  final  quantity  of  interest  is  given  by  (5)  in  terms  of  |C|4.  For  the  sum  in  (6),  the  latter  quantity  is 
given  by 

fc!  ~  k^~  **  *■  (13) 

This  statistic  is  evaluated  in  appendix  A.  In  fact,  the  more  general  quantity 

b|  CM  (14) 

is  evaluated,  where  complex  random  variables  ak,  bk,  ck,  dk  are  statistically  dependent  amongst 
themselves  for  any  k;  but  these  random  variables  are  statistically  independent  of  the  random  variables 
an,  bn,  cn,  dn  for  all  n  *  k.  We  evaluate  (14),  rather  than  (13),  for  two  reasons:  First,  the  notation  is 
simpler  and  an  error  in  analysis  is  much  easier  to  detect;  secondly,  this  more  general  result  may  have  a 
possible  future  application,  and  the  analytical  effort  is  no  greater.  The  average  in  (13)  is  given  by  (A- 
26)  in  terms  of  the  16  fundamental  sums  defined  in  (A-24),  and  the  variance  of  |C|2  is  given  by  (A-31). 
We  do  not  repeat  these  results  here  because  of  their  length.  A  program  for  evaluating  all  the  above 
quantities  is  given  in  appendix  A,  table  A-l . 

As  a  check  on  these  results,  the  case  of  Gaussian  complex  random  variables  {zk}  is  considered  in 
appendix  B.  A  program  for  this  special  case  is  given  in  appendix  B,  table  B-l . 

APPLICATION  TO  BEAMFORMING 


For  an  ideal  array  with  no  perturbations  in  element  positions,  gains,  phases,  delays,  or  failures,  the 
voltage  transfer  function  to  a  plane-wave  arrival  can  be  expressed  as 


where  vk  is  a  complex  quantity  that  incorporates  the  plane-wave  arrival  angle  and  frequency,  the 
steering  angle,  and  the  element  parameters  such  as  element  position,  gain,  and  phase;  see  reference  1, 
equations  (1)  -  (8).  When  imperfections  in  the  array  realization  are  encountered,  they  can  be  included 
in  the  array  voltage  transfer  function  by  replacing  vk  in  (15)  by 


(16) 


where  gk,  rk,  $k  are  real  random  variables.  The  random  variable  gk  represents  random  element 
failures,  by  setting  gk  =  0  or  1  with  specified  probabilities;  the  random  variable  rk  represents  relative 
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gain  perturbations  from  the  desired  value  of  zk  =  vk;  and  random  variable  $k  represents  phase 
perturbations  from  the  desired  value  of  0. 

The  phase  perturbations  {<j>k}  can  arise  from  positional  and/or  delay  perturbations  in  the  array 
realization.  This  problem  is  considered  in  appendix  C  for  a  linear  array;  the  variance  of  4>k  is  derived 
and  its  dependence  on  arrival  and  look  angles  is  made  explicit  for  two  situations  of  knowledge  on  the 
part  of  the  array  designer. 

The  random  variables  gk,  rk,  $k  are  dimensionless.  We  shall  assume  that  they  are  independent  of 
each  other  and  of  all  other  random  variables  for  different  values  of  k.  (This  could  be  generalized  at 
the  expense  of  requiring  more  detailed  knowledge  of  the  joint  statistics  of  these  random  variables.) 
We  let  the  moments  of  these  random  variables  be  denoted  by 

J  M) 

of  k . 

E  X.  (<®*?4  _  (17) 

The  independence  of  these  moments  on  k  could  be  generalized  easily,  but  is  not  done  herein.  Thus, 
physically,  we  are  presuming  an  array  where  all  elements  are  equally  random  in  terms  of  amplitude 
perturbations,  phase  perturbations,  and  failures. 

The  necessary  statistics  that  must  be  evaluated  are  listed  in  (A-25).  They  are,  using  (16)  and  (17), 
given  by 

%  » 

K  =  VK  ^ 

W  (18) 


It  is  worthwhile  noting  that  only  the  first  two  moments  and  y2  are  required  for  the  phase  per 
turbation.  Then  the  fundamental  sums  that  {Ty}  in  (A-24)  depend  on  are  just 


which  are  independent  of  the  perturbation  statistics  in  (17).  The  reason  for  this  independence  is  that 
(17)  was  presumed  independent  of  k,  the  element  number. 


Now  we  can  express  the  various  array  responses  of  interest  in  terms  of  the  above  quantities.  We 
have  from  (15)  (or  from  (16)  for  gk  =  l,rk  =  0,  ^k  =  0,  all  k) 


Ideal  (Complex)  Voltage  Response  =  ^  v. 


(20) 
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Then 

Ideal  Power  Response  =  |^vk|  .  (21) 

Next,  from  (16), 

Actual  Voltage  Response  =  C  *  .512,  (22) 

and  finally 

Actual  Power  Response  =  )C  I* .  (23) 

Then  (22),  (7),  and  (18)  yield 

Average  Voltage  Response  =  Av[0 1  =  ,  (24) 


which  is  a  scaled  version  of  the  ideal  voltage  response  (20).  This  simplification  results  because  the 
moments  in  (17)  were  assumed  independent  of  element  number  k.  From  (23),  (9),  (18),  and  (24),  there 
follows 

Average  Power  Response  =  Av  l  )C|*} 

•  (25) 

the  second  term  of  which  is  a  scaled  version  of  the  ideal  power  response  (21).  The  first  term  of  (25)  is 

the  variance  of  the  voltage  response,  as  may  be  seen  by  combining  (22),  (10),  and  (18): 

Variance  of  Voltage  Response  =  Var  {  c) 

=  (/»**-;*>,*  Ilf,  (26) 

Finally,  the  quantity 

Variance  of  Power  Response  =  Var  [|  C  (27) 

is  given  by  (A-31),  (A-24),  and  (18). 

EXAMPLE  OF  ELEMENT  PERTURBATIONS 

The  necessary  moments  were  listed  in  (17).  We  now  need  to  specify  the  probability  density  func¬ 
tions  of  (gk),  {rk},  and  {$k}  in  order  to  evaluate  jim»  vm,  rnv  Since  gk  is  a  0,  1  random  variable 
representing  element  failures,  we  let 

p(g)  -  Q  S(9>  +*  0-  Q)  V  ($'•)  i  (28) 

that  is,  Q  is  the  probability  of  element  failure.  We  have  taken  advantage  of  the  independence  of  k  in 
(17).  Then 

fr*”  Y  9  f  all  (29) 
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The  second  quantity  in  (17)  is 

)>«  =  0  +  r)U  (30) 

If  we  let  r™  *  ftm,  then  we  have 

»,  -  *4  *  ^  *  I+Y+Y+yy*-  •  (31) 

As  a  special  case,  if  relative  perturbation  r  is  zero-mean  Gaussian,  then 

fW-~re»r(-^r-)  .  02) 

and  there  follows 

p,  =  0,  ^  *  <r;,  =  o,  ^  =  3<rr\  (33) 

for  which 


>>,  =  l,  4  -  1+  <t;\  ■  1+  3(rra  ,  =  |+  6  rr*  +  3r* .  (34) 

Finally,  to  evaluate  the  third  quantity  in  (17),  we  need 

^  ■  ;#.*>#).  ,35» 

The  probability  density  function  of  $  will  be  taken  to  be  zero-mean  Gaussian,  in  which  case, 

X,  *  J 4*  e*f  (- f£)  *  ***  (-  >  m*).  (36) 

It  should  be  noticed  that  ym  is  real  for  this  example;  this  real  property  simplifies  the  programming 
eff  rt  and  is  used  throughout  the  rest  of  this  report. 


EQUISPACED  LINEAR  ARRAY 


GENERAL  RESULTS 

It  was  noted  earlier  that  the  voltage  transfer  function  of  an  ideal  array  to  a  plane-wave  arrival  can 
be  written  in  the  form  of  (15).  We  now  investigate  this  form  for  an  equispaced  linear  array;  we  find, 
for  a  symmetric  real  weight  structure  { wk }  about  the  center  of  the  array,  that  (reference  1 ,  page  3) 

VK  *  V*n  i  (k-  4)u"3  ’ffr  on  ever*  Pim«W  of  elements,  (37) 


where 


U  -  2-n{  -f  (  Sin  ^  -  Sin  =  Zir  ^  -  5i»  ^  (38) 

Here  fa,  c,  A,,  and  arc  the  frequency,  speed  of  propagation,  wavelength,  and  arrival  angle 
(measured  from  broadside)  of  the  plane-wave  arrival;  d  is  the  linear  array  element  spacing;  and  is 
the  look  (steering)  angle  of  the  array.  Extension  of  (37)  to  an  odd  number  of  elements  is  readily 
achieved,  but  not  pursued  here.  We  let  N  be  the  total  number  of  elements  in  the  array  and  express 


6 
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N  -  2H,  (39) 

where  H  is  the  number  of  elements  in  one-half  of  the  array. 

The  five  fundamental  sums  in  (19)  now  take  the  form 

coa[(k-i)«]  «  L,(«) 

fKf-  V, 

=  2L  <  as[(2IH)a]  «  L,  (u) 

^Kl\  •  cbs[(|<-4)u]  *  L» 

4^"  *  *TV5,  (40) 

where  we  have  taken  advantage  of  the  symmetric  real  weight  structure  in  order  to  write  all  the  sums  as 
explicitly  real  quantities.  The  notation  I  here  denotes  a  sum  over  all  nonzero  weights  {wk}  from 
k  *  -  H  to  H.  In  the  program  to  be  presented  later,  advantage  is  taken  of  the  symmetry  in  order  to 
decrease  the  number  of  terms  computed  by  a  factor  of  2.  If  we  denote  the  first  function  in  (40)  by 
L,(u),  then  it  follows  that 

L,(-m)  -  L|  («),  L,  (tt+2rr)  *  -  L,  («).  (41) 

Similar  useful  properties  hold  for  the  other  L  functions  in  (40);  they  enable  the  region,  where  (27) 
must  be  computed,  to  be  reduced  to  the  range  (0,n). 

In  terms  of  the  quantities  defined  in  (40),  we  can  now  express  (20)-(26)  as 

X  VoHtpji  KeS ponj»  ■■  L,  (u) 

Xdfa)  Po vy<tr  «  L*(u) 

A vwpfle  VeHoje  } *  C,  L,  (u) 

Avcynje.  fta*v  J?efj>onje  •  C*  Lf (t«)  +  (C**- 

Vdtiawfv  *f  Vo) tnjt  Rejpowse  *  ,  (42) 

where 

,  C2„  *  ■  (43) 

The  variance  of  the  power  response  is  given  in  appendix  D.  Examples  of  (42)  and  the  variance  of  the 
power  response  are  deferred  until  later  in  this  report.  A  program  for  calculating  the  average  behavior 
for  a  linear  array  is  given  in  appendix  D,  table  D-l , 
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DEEP  SIDELOBE  BEHAVIOR 


In  the  deep  sidelcbe  region,  the  ideal  voltage  response  is  substantially  zero;  that  is,  L,  (u)  in  (42)  is 
approximately  zero.  Then  (42)  simplifies  to 

Vo)fc»y*  'fcjpOl*#  •  0 

Xdw/  Poiotr  R«yp  t>*Dt  *,  o 
7frrpo«e  *  o 

A>er«}e  potv«r  p4  ~  jif  y*)  iw * 

Var )'«»«  of  Vo  H»Jt  Potpww*  «  ()4A  -  /tf  Wf)  %-*£  ,  (44) 

where  we  employed  (43)  and  (40).  The  variance  of  the  power  response  is  given  in  (D-6)  in  terms  of  the 
other  quantities  defined  there;  no  simple  expression  for  this  variance  was  obtained.  Alternative  in¬ 
terpretations  of  the  constants  are  given  in  appendix  D,  (D-9). 


Since  the  peak  ideal  power  response  (for  positive  weights)  is  given  by  (see  (42)  and  (40)) 

L*  (0)  *  (frf  . 

then  it  follows  from  (44)  and  (45)  that 

AVe*a<f  R >m«r  foapows*  (dop  region)  =  m  V_ 

Ptflfc  I<i«o!  Power  P%»pow»e  Ntff 


where 


N 


iff 


M. 


^w* 


(45) 


(46) 


(47) 


The  quantity  in  (46)  requires  only  second-order  statistics  of  the  perturbations  and  the  single  summary 
parameter  (47)  of  the  weights.  Equation  (47)  is  maximized  by  equal  weights  over  the  entire  array. 


For  the  example  considered  earlier,  we  use  (29),  (34),  and  (36)  to  evaluate  the  numerator  of  (46)  as 

V  ■  0-Q)(l+  *r)  ~  (l  “  <?)*  e*f  (-  <).  (48) 

For  small  phase  perturbations,  a2p«  1,  and  (48)  simplifies  to 

V  *  (M)  (9  4  r*  ♦  ff'/O-Q)).  (49) 

If,  additionally,  the  probability  of  element  failure  is  small,  Q  «  1 ,  then  (49)  can  be  further 
manipulatedinto  a  variety  of  forms: 

V  *  (l -Q)  (<?■*■  r’  +  'rr‘) 

*  0 -<*)$-*-( 

*  <T^  +-  ifr  f  Trek (eW»wt  okoj) 

s  r*  +  r?  +  <  ■  (50) 


4 


8 


TR6191 


Then  (46)  becomes 

fW  ftnww (Ao  «sj»lyk  .  /<i\ 

Jti),  XckJ  Rw*r  Ki tf*»  N,^ 

This  simple  rule-of-thumb  is  similar  to  (43)  in  reference  1  and  to  (83)  in  reference  2,  but  now 
generalized  to  include  simultaneous  element  failures  and  gain  and  phase  perturbations.  It  indicates 
that  the  variances  of  the  three  perturbations  simply  add,  at  least  under  the  assumptions  made. 

GAUSSIAN  APPROXIMATION  IN  DEEP  SIDELOBE  REGION 

When  the  number  of  elements  in  the  array  is  large,*  the  sum  variable  C  may  be  well  approximated 
as  a  complex  Gaussian  random  variable.  However,  the  statistical  description  of  C  is  still  difficult, 
since  the  variances  of  the  real  and  imaginary  parts  of  C  are  not  equal,  in  general;  see  reference  1, 
appendixes  A  and  B.  However,  in  the  deep  sidelobe  region,  we  have  from  (22),  (24),  (25),  and  (40), 

C  *0,  Cr*  O 

lef  #  V) •  (52) 

Then,  letting  complex  random  variable 

(53) 

we  have 

~x  -  y  «  o ,  xy  *  0,  xT  -  Y  ■  »  •  (54) 

Now,  if  none  of  the  elements  have  failed,  the  deep  sidelobe  behavior  is  attainable,  and  we  have, 
under  this  conditional  situation,  the  joint  Gaussian  probability  density  function 

e*p(“  >  <55> 

where,  from  (52), 

<  =  ,  (56) 

since  we  must  set 

=  i  ftr  all  m  (57) 

under  this  condition  of  no  elements  failing.  Then 

As  [Icfl  - 

Vo r  f I  e|^  *  Vo rix%</}  *  r? 

3W  Pev  {|C|^  =  r*.  (58) 

•The  material  in  this  section  applies  to  any  array,  not  just  a  linear  array;  in  fact,  the  larger  the  number 
of  elements,  the  better  the  approximation. 
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Thus,  the  standard  deviation  of  power  gain,  |C|Z,  equals  its  mean;  this  property  obtains  only  when 
none  of  the  elements  have  failed.  Also 

ftok  (|C|’<  T)  -  FU  (s* <  t)  *  I  -  t'xj (-T /<rt‘) .  (59) 

That  is,  the  probability  density  function  of  |C|2  is  exponential  in  the  deep  sidelobe  region.  In  par¬ 
ticular, 


Avfc|']+k*5{4&v{lci1)  «  I  -  %*}>(- i-k).  (60) 

Thus,  for  example,  the  probability  that  the  actual  power  gain,  |C|2,  remains  less  than  its  average  value 
is  0.63,  while  the  probability  that  it  is  less  than  the  mean  plus  two  standard  deviations  is  0.93. 

PLANAR  ARRAY 


GENERAL  RESULTS 

In  order  to  apply  the  results  above  to  a  planar  array,  it  is  convenient  to  use  a  one-to-one  number 
association  between  the  integer  k  and  the  location  p,  q  of  a  particular  element  in  a  planar  array 
presumed  to  have  a  grid  structure.  Thus,  instead  of  (16),  we  have 


(61) 


and  instead  of  (17),  we  have 


■  A 


hrJ  *  a 


itldepe*d»i^'  of-  p,  y 


J 


(62) 


The  necessary  statistics  that  must  be  evaluated  are  identical  to  (18),  except  that  vk  is  replaced 
everywhere  by  \pq.  Thus,  analogous  to  (19),  we  need  to  evaluate  the  quantities 

W  »  ^  vn  >  • 


, 


n  "  ft 

Analogous  to  (20)  et  seq.,  we  have 


Xc)e«l  (Co^W)  Vo  Response  *  ^  Vft 


for  gpq  ■  1.  fpq  -  o,  ^  -  0,  all  p,q.  Then, 


Next,  from  (61), 


Xdtoi  Vomer  Vespome  “  | 


Ackw\  Voyage  F?espo«J«  -  Cr  =  2 

Tv 


(63) 


(64) 


(65) 


(66) 
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jAc”tu«|  Ru*r  ■  |CJ- .  (67) 

Then,  (66),  (7),  (A-24),  and  (62)  yield 

A )?ejpon»*  •  Av[c)-  t,  *  /.“.V,  ,  (68) 

which  i3  a  scaled  version  of  the  ideal  voltage  response  (64).  This  simplification  results  because  the 
moments  in  (62)  were  assumed  independent  of  element  location  p,  q.  Also,  from  (67),  (9),  (A-30),  (A- 
24),  (61),  and  (62),  there  follows 

Avev^je  Powwr  Re3j»o*,«  -  Av  [|C|*}  »  T*,-  "Q,  +•  }x|m 

.  +  )  ^n\\  m 

the  second  term  of  which  is  a  scaled  version  of  the  ideal  power  response  (65).  The  first  term  of  (69)  is 
the  variance  of  the  voltage  response,  as  may  be  seen  by  combining  (66),  (10),  (A-24),  (61),  and  (62): 

Mari  once  of  V«K«fo  Kotfows*  *  \/or  \  c]  * 

*  M.  (70) 


Finally  the  quantity 


Vovianr«  of  Pohmt  Reifo**  ■  Vor^lCl*} 


is  given  by  (A-31),  (A-24),  and  (18)  with  vk  replaced  by  vN  everywhere;  the  quantities  that  must  be 
evaluated  are  those  listed  in  (63). 

EQUISPACED  PLANAR  ARRAY  WITH  MULTIPLICATIVE  WEIGHT  STRUCTURE 

For  a  planar  array  with  elements  equispaced  on  the  x,  y  plane  by  distances  dx,  dy,  and  for  a 
multiplicative  weight  structure, 


w  -w“w»  , 


we  have,  from  (21)  and  (22)  of  reference  1 , 


% "  w? e*F  [-  •  (p  *  ^ 3  «?[-•  (r 


for  a  planar  array  with  an  even  number  of  elements  in  both  the  x  and  y  coordinates,  where  the  weight 
structures  {w£x>}  and  {w^>}  are  assumed  symmetric  about  both  x  *  0  and  y  =  0  (p  =  q  =  0),  the 
center  of  the  array.  The  parameters  u  and  v  incorporate  look  (steering)  angle  ( 9 ,,  f,),  spacings  d*,  dy, 
and  plane-wave  arrival  wavelength  A,  and  angle  (0t,  $a): 
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u  »  2ir-^  (a'n  4  ^  ~  5Wf» t<>3  ©i) , 

V  *  2tr  (sin  ^  «'»*<$-  si«4,  *'"  8,), 


(74) 


The  polar  angle  (measured  from  the  z-axis)  is  4,  and  the  azimuthal  angle  (measured  from  the  x-axis)  is 
6.  The  dimensionless  parameters  required  in  (74)  are  relative  spacings  d/A,,  d/A,,  arrival  angle  (0a, 
4a),  and  look  angle  (0(,  4i)< 

The  five  fundamental  summations  required  in  (63)  take  the  form 
Zvn-  c.s[(>-  *)*]•  4)v] 

2W-  *xww? 

C03 [(2?-t)u]  •SwJ)'  fe>s[fy- 1)  v)  t  L?(4L?(v) 
r>  r  V 

2|vn|’V  X'v,*’ r<»[(p-d^  c®s[(^-  J[)v]  s  L*J  (w)  L^C») 

»r  f  4 

Slvr  ■  Z.<-  2^  *vN!*. 

n  n  p  f  ^  (75) 

These  quantities  are  all  real,  due  to  the  symmetry  assumptions;  so,  for  the  perturbation  example  in 
(28)-(36),  all  the  {T^}  in  (A-24)  are  real.  Furthermore,  the  summations  on  negative  p  and  q  in  (75)  can 
be  avoided  by  multiplying  the  positive  p  and  q  terms  by  a  factor  of  2,  as  done  earlier  for  the  linear 
array. 

We  can  express  the  desired  quantitites  in  (64)-(71)  in  terms  of  (75)  and  (43): 

X4««l  Reip  0«fC  *  Lpw  l“w 

Idr«)  PwHtr  Pejponje  *  L^(y)  (v) 

Avojra^e  Vo  )to0e  Rei^o^ie  -  C,  L^(u) 

frvroqt  po«*ur  Rejpome  «  (c*.- C*)W’k^+  Cf  L^llf^v) 

Vw  i  a*«  of  VoKage  Pespo^je  *  C*)'Wf>'W^>.  (76) 

The  variance  of  the  power  response  is  given  by  (D-l  1)  in  appendix  D.  A  program  for  calculating  the 
average  behaviors  for  the  planar  array  is  given  in  appendix  D,  table  D-2. 

EXAMPLES 

Four  curves  are  drawn  in  each  of  the  figures  discussed  next.  The  bottommost  curve  in  the  deep 
sidelobe  region  (the  curve  with  the  deep  notches)  is  the  ideal  power  response  (21),  which  would  be 
realized  for  no  element  failures,  gain  perturbations,  or  phase  perturbations;  this  curve  is  normalized 
to  0  dB  at  its  peak  where  the  look  angle  equals  the  arrival  angle.  The  second  curve  from  the  bottom  (in 
the  deep  sidelobe  region)  is  the  average  power  response  (25)  for  the  particular  set  of  perturbation 
statistics  listed  with  each  figure.  The  third  curve  from  the  bottom  is  a  plot  of  the  average  plus  one 
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standard  deviation  of  the  power  response,  that  is,  (25)  plus  the  square  root  of  (27).  The  topmost  curve 
is  the  average  plus  two  standard  deviations  of  the  power  response,  that  is,  (25)  plus  two  times  the 
square  root  of  (27).  In  the  mainlobe  region,  the  curves  can  cross  each  other;  however,  the  ideal  power 
response  curve  always  reaches  0  dB,  and  the  other  three  average  curves  always  lie  in  the  order  relation 
indicated  above.  These  facts  enable  the  reader  to  discern  the  behavior  of  the  four  curves  in  any  one 
figure. 

LINEAR  ARRAY 

We  consider  an  equispaced  linear  array  of  N  ■*  20  elements.  The  single  variable,  that  all  the  array 
responses  depend  on,  is  the  variable  u  defined  in  (38), 

a  *  2^^  (sin  .  <77> 

for  which  the  range  (0,n)  is  sufficient  to  cover  ail  cases  of  element  spacing,  look  angle,  and  arrival 
angle  and  wavelength. 

In  figure  1,  a  Dolph-Chebyshev  -30  dB  array  design  is  indicated,  for  which  the  effective  number  of 
elements,  (47),  is  equal  to  17.35.  This  result  is  especially  useful  for  evaluating  the  average  power 
response  in  the  deep  sidelobe  region;  see  (46)  and  (51). 

The  four  parts  (A,  B,  C,  D)  of  figure  1  correspond,  respectively,  to  (A)  element  failures  only;  (B) 
weight  perturbations  only;  (C)  phase  perturbations  only;  and  (D)  combined  element  failures,  weight 
perturbations,  and  phase  perturbations.  Quantitatively,  we  have 

(A)  probability  of  element  failure  Q  •  0.001 

(B)  variance  of  relative  weight  perturbation  oj  m  0.001 

(C)  variance  of  phase  perturbation  o£  =  0.001 

(D)  all  the  above  combined . 

The  parameter  values  have  been  chosen  so  that  the  three  variances  are  equal 
(Og  "  Q(1  ~  Q)  a  0.001);  thus  the  average  power  response  in  the  deep  sidelobe  region,  (51),  should  be 
equal  for  parts  (A),  (B),  and  (C).  In  fact,  (51)  gives  0.001/17.35  =  -37.6  dB  for  the  first  three  parts 
of  figure  1  and  3  x  0.001/17.35  -  -42.4  dB  for  part  (D).  These  calculations  agree  very  well  with  the 
results  plotted  in  figure  1. 

The  curves  in  figures  IB  and  1C,  for  weight  and  phase  perturbations,  are  virtually  identical;  the 
curves  in  figure  1  A,  for  element  failure,  indicate  slightly  poorer  performance,  about  1  dB  larger  at  the 
peaks  of  the  sidelobes.  Figure  ID,  for  combined  perturbations,  is,  of  course,  the  poorest  of  all.  The 
mainlobe  response  is  substantially  unchanged  in  the  four  parts  of  figure  1 . 

Figure  2  is  drawn  under  conditions  identical  to  figure  1  except  that  the  weighting  is  changed  to 
Hamming.  For  this  N  *  20  element  array,  Neff  -  14.68.  The  results  for  the  average  power  response  in 
figures  2A,  2B,  and  2C  are  virtually  identical;  however  the  standard  deviation  of  the  power  response 
for  element  failures  is  somewhat  greater,  thereby  leading  to  poorer  performance  in  figure  2A. 

Figure  3  is  also  drawn  under  identical  conditions  except  that  the  weighting  is  now  Hanning,  for 
which  Npff  *  13.33.  The  observations  for  this  figure  are  identical  to  those  for  figure  2. 
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Figure  1A.  Q  -  .001,  oj  •  0,  o|  ■»  0 


Figure  IB.  Q  -  0,  o*  -  .001,  o*  -  0 


Figure 1.  Equlspaced  Linear  Array;  -30  dB  Dolph-Chebythev  Weights 


Figaro  ID.  Q  -  .001,  oj  ■  .001,  oj  ■  .001 


Figure  1.  (Coat'd)  Eqalspacod  Llaoor  Array;  -30  dB  Dolpb-Chobyslrev  Weights 


Figure  3A.  Q  «  .001,  o?  “  0,  cj  «  0 


Figure  38.  Q  -  0,  •  .001,  oj  ■  0 


Figure  3.  Equispuced  Linear  Array;  Hanning  Weight* 
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Figure  3C.  Q  »  0,  of  *  0,  of  «  .001 


Figuc  3D.  Q  -  .001,  of  -  .001,  of  -  .001 

flfur*  3.  (Coat’d)  EqoJspacad  Uaw  Arvay;  Haaaiaf  Weight* 
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PLANAR  ARRAY 

We  consider  an  equispaced  planar  array  of  Nx  x  Ny  *  20  x  18  elements,  for  a  total  of  360 
elements.  Since  there  is  no  simple  single  parameter  like  u,  (77),  in  the  linear  array,  we  plot  the  results 
versus  the  polar  arrival  angle  or  azimuthal  arrival  angle,  for  specified  polar  and  azimuthal  look 
angles.  The  results  in  figure  4  et  seq.  employ  element  spacings  in  the  x  and  y  coordinates  of  a  half 
wavelength;  however  the  program  in  table  D-2  easily  accommodates  other  spacings  of  interest.  Figure 
4  is  drawn  for  a  multiplicative  weighting  of  -30  dB  Dolph-Chebyshev  design  in  each  of  the  x  and  y 
coordinates;  the  effective  numbers  of  elements  in  the  x  and  y  coordinates  are  17.35,  i5.57,  respec¬ 
tively.  The  average  responses  in  figure  4  (notice  the  much  larger  element  perturbations)  are  plotted 
versus  the  polar  arrival  angle  $a  in  radians,  while  the  azimuthal  arrival  angle  and  the  look  angles  are 
all  zero.  The  results  in  figures  4A,  4B,  and  4C  are  substantially  the  same;  the  case  for  combined 
perturbations  in  figure  4D  is,  of  course,  poorer  in  terms  of  performance  attainable. 

The  average  power  response  at  the  peak  of  the  mainlobe  has  dropped  by  about  1  dB  in  figure  4A 
because  of  the  effect  of  element  failures,  and  by  about  1/2  dB  in  figure  4C  because  of  the  effect  of 
phase  perturbations.  The  combined  effect  in  figure  4D  is  such  as  to  lower  the  average  power  response 
at  the  peak  by  about  2  dB. 

The  similarity  of  results,  for  weight  perturbations  or  phase  perturbations  or  element  failures  alone, 
prompts  us  to  confine  further  attention  only  to  the  case  of  combined  nonzero  perturbations.  In  figure 
5,  the  same  array  is  used  as  in  figure  4,  where  now  we  have  selected  Q  =  o*  *  o*  *  0.1.  Part  A 
corresponds  to  a  plot  of  average  results  versus  azimuthal  arrival  angle  0a,  while  part  B  is  for  a  varying 
polar  arrival  angle  $a.  The  reason  that  the  ideal  power  response  goes  below  -60  dB  is  a  result  of  the 
multiplicative  effect  of  the  individual  x  and  y  -30  dB  patterns.  This  leads  to  average  responses  that  are 
virtually  independent  of  the  polar  and  azimuthal  arrival  angles. 

Figure  6  is  drawn  under  conditions  similar  to  figure  5  except  that  the  weighting  is  changed  to 
Hanning.  The  effective  numbers  of  the  20  x  18  elements  are  13.33  x  12  in  the  x  and  y  coordinates. 
Finally,  in  figure  7,  when  we  make  the  polar  look  and  arrival  angles  <^,  4>a  equal  to  n/2,  we  get  two 
large  equal  responses  at  0a  =  0,  and  +  n.  This  is  a  result  of  the  element  spacings  being  equal  to  a 
half  wavelength  and,  therefore,  the  array  is  unable  to  distinguish  or  reject  arrivals  coming  endfire  to 
the  array  from  the  undesired  opposite  direction. 

Additional  results  for  other  equispaced  linear  or  planar  arrays  are  easily  available  by  using  the 
programs  in  appendix  D  for  whatever  set  of  parameters  fits  the  user’s  application.  Generalizations 
should  be  obvious  and  easily  realized  from  the  general  results  presented  in  appendix  A. 

SUMMARY 

General  equations  for  the  first  four  moments  of  a  sum  of  independent  complex  random  variables 
have  been  derived  with  no  restrictions  on  the  statistics  of  the  individual  variables  and  no  requirement 
of  identical  statistics.  A  program  for  this  case  is  available  in  appendix  A,  table  A-l,  for  those  ap¬ 
plications  requiring  the  most  general  case.  This  program  has  been  thoroughly  checked  and  also 
compared  with  the  results  for  Gaussian  random  ’'ariables. 

These  general  equations  have  been  specialized  here  to  a  beamforming  application,  including 
equispaced  linear  and  planar  arrays.  Additionally,  for  ease  of  programming  and  investigation, 
identical  statistics  were  assumed  for  all  the  array  elements;  however  this  restriction  could  easily  be 
eliminated  by  reference  to  the  general  results  given  earlier.  To  determine  allowable  tolerances  on 
element  weight  perturbations,  phase  perturbations,  and  element  failures,  it  is  recommended  that 
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initial  guesses  be  made  at  the  statistical  parameters  and,  then,  the  average  responses  given  herein 
plotted  and  observed.  Then  another  educated  guess  can  be  made  and  the  new  responses  observed.  In 
this  manner  one  can  quickly  converge  to  tolerable  limits  on  the  different  types  of  perturbations  in 
order  to  realize  specified  performance  and  sidelobe  levels. 


FlgeredA.  Q  «  .1,  o*  «■  0,  oj  »  0 


♦. 


Figure  4B.  Q  »  0,  oj  »  .1,  oj  *  tf 

Figure  4.  Equbpeced  Pfaunr  Array;  -30  <1B  Dolpb-Cbebyshev  Weights,  ^  -  0,9t  ■  0 
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Figure  4C.  Q  «  0,  o*  •  0,  a*  -  .1 


Figure  4D.  Q  »  .1,  of  “  .1,  oj  *  .1 


Figure  4.  (Coat’d)  Eq  ah  paced  Planar  Array;  -30  dB  Doiph*Cbebyshe?  Weights,  +j  -  0, 0,  •  0 
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Fl*mi«  7.  Eqaiapaccd  Pl&sar  Array,  Hutting  Weights,  -  n/2,  Ot  »  n/t,  f,  »  r/2,  Q  » 
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,1,  of  •  .1,  o|  *■  .1 
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APPENDIX  A 

DERIVATION  OF  FOURTH-ORDER  MOMENT 


GENERAL  RESULTS 
Here  we  will  evaluate  the  complex  quantity 

1  (A-» 

where  complex  random  variables  ak,  bk,  ck,  dk  are  statistically  dependent  amongst  themselves  for  any 
k,  but  these  random  variables  are  statistically  independent  of  the  random  variables  an,  b„,  cn,  dn  for 
all  n  *  k.  The  fourth-order  sum  in  (A-l)  can  be  broken  down  into  many  subcases:  (1)  all  four  sub¬ 
scripts  equal;  (2)  three  equal  and  one  different;  (3)  two  pairs  of  two  equal,  but  not  to  each  other;  (4) 
two  equal  and  two  others  different  from  each  other  and  the  two  equal  ones;  and  (3)  all  different.  The 


values  of  the  components  of  S  in  each  subcase  are 

S  *  S.+-  5af  V  5*  +  5r  ,  (A-2) 

where 

5,  *  ,  (A-3) 

(A-4) 

S3  •  £  4  ^ir  4  !>».(■,  }  (A-5) 

S+  -  c* +  £ a*  ^  ^  5  c« 

4  4  4  >  (A-6) 

s,  =  X  0*  5  s  £  »  (A-7) 


and  where  T  denotes  summation  over  all  subscripts  that  are  unequal  in  the  summands.  (As  a  check, 
the  number  of  terms  in  S,-S5  are  N,  4N(N-1),  3N(N-1),  6N(N-l)(N-2),  and  N(N-l)(N-2)(N-3), 
respectively;  these  add  to  N4,  as  they  must.) 

The  summations  in  (A-3 HA-7)  must  be  simplified.  To  do  so,  we  adopt  an  abbreviated  notation; 
for  example,  (A-3)  is  denoted  as 


5,  =  *  JBCET 

K 


(A-8) 


A-l 
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In  order  to  simplify  (A-4)  and  (A-5),  we  develop 


f-v4  ■«  UV-UV. 

(A-9) 

Notice  that  U-V  denotes  the  product  of  two  summations,  whereas  UV  denotes  one  summation.  Then 
(A-4)  becomes 

+  *&B->™+15CB*X-FcdA, 

(A-10) 

while  (A-5)  yields 

S3  •  AS-  c»  -  AkB  +  Acbd-acbd  +  Ad-bc-ABbc  . 

(A-U) 

Thus,  for  example,  the  first  two  terms  in  (A- 10)  are 

AK  D -A&cD  .  ■ 

(A-12) 

To  simplify  (A-6),  we  use  the  development 


^  Sx  •  J 

V  K  n^k  tit*  k  ««»  L**  J 

--  W0>vC-  -  5-u.kv. 

j  o  u 

=  w-(uv-uv)~  (u-vw-uvw)-  (uwv-uvw) 

:  VJ  V-W-  IJV-W-  1>VW-  V-Utf  +  2  UVW,  (A-13) 

where  we  used  (A-9).  Then  A-6  becomes 

v+  --  AB  OD-ABc-E- AB  CD-  C-AfP  +  2A105 
+  AC  B  D-ACB  D  -  Ac-  Bp  -  &AC0  4-  2AZ65 
+  AP-B-C- At>B-C  Ad-BC-B-AK  *7APBC 
+  bc-a-B-Abc  5  -Bc  Ap-A-Bc^  +  2BcA& 

+  BC-A  c-  Abp-C-BC-Ac-  A-bBc  +  20£Ac 

4  tp.A.E-Afp-B-  C6-M-AF-BC5  +  2CBAB.  (A-14) 

Finally,  we  need  the  development 


i 


! 


i 


( 


! 

i 


A-2 
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?“****  ’  |"S#*  §  t^VW^] 

=  X'£i\W  Z.Un’iV.K 

*  X'  [u  v  w-  UV  W-  U-VW  *  V  WW  *  2  uvw] 

-l/X-V -W+  UVX-W4  UX*VW  +V  UWX-  2UVWX 
-VX'U  W  +  UVX-W4-  UVWX  +  XV-UW-  2UVWX 
-  Vvx  o-V  +UV  WX+ U'VWX  +  V'UWX-iUVwy 
«  U-VW-  x~  (w  XW^-VW'UX+UWX- V+  UX-  VIV  4 VX‘UW+WX'UV) 

+  2  (uvw  X  +  uwxv  +  uvx  w  +  vwx-  u) 

4-UX  VW+  XV  UW  +  UVWX-  (  UVWX,  (A-15) 

where  we  used  (A-13).  Then  (A-7)  becomes 

^  -A'g  CS-  (tec  D+ICA-D+AO  §D4  AD  8-c  +  BM-  c  +  cd-*B) 

4- 2(A-6c  6 4 B  +  >^BD  C  +  Bc6 

+  Ab-c6  4-  #C-  65  +  *&•  6  c  -  Ubcd  .  (A- 1 6) 

It  is  not  possible  to  simplify  or  combine  any  of  the  60  terms  in  (A-8),  (A-10),  (A-ll),  (A-14),  and 
(A-16)  because  they  each  employ  different  statistics.  Thus,  the  general  answer  to  (A-l)  is  given  by  (A- 
2),  where  the  five  components  are  given  in  the  equations  just  listed  above.  Repeating  the  notation,  the 
first  and  last  terms  in  (A-16)  are  given  by  ak)  Bf)  (Z  cm)  (Z  3n)  and  -6  ^  ak  6k  ck  3k,  respectively. 

SPECIAL  CASE 
We  now  set 

(A-l  7) 

in  which  case  (A-l)  specializes  to 

S  «  2-  _  1 }  (A-l 8) 

that  is,  S  is  now  the  mean  value  of  the  magnitude-fourth  power  of  the  sum.  Then  (A-8)  yields 

s,  -  TzT  -  ?  kp  •  (*-'»> 

(A-10)  yields 


A-3 


'T’*'**  ItWW'  TfWlwWB  *’**  *  ”'WT* 


”  ^TOTSi^Wqn^  ■* 
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5, » 2|?P?  '2*  -21177  +  comr^  <*"vM3*Jk 
.  4Rc(izf2'f 

-fwstf]. 

and  (A-i  1)  becomes 

s,.  jiiT  iiP-^iJi*  +  ?■  ?-  1ST 

.2(fB*f-2|B#  +  |?*r-  ?ISI- 


(A-20) 


(A-21) 


The  quantity  S4  in  (A-14)  can  be  developed  as  follows: 

s, .  ♦{  w-  z- i-  lip  ??- TzT-  i?r-w?’?  +^6Fi*rj 

42ll{  +  2Zr(fJ} 

2fc{f E?S  Sf)  -  f>5F-2KI‘*25iir lil'] 

+  2 -225 +  (A-22) 

Finally  (A- 16)  yields 

tisr 

4  2x}+2(|.isi^4|^|'-  <?  W\  (A.23) 

At  this  point,  it  is  convenient  to  define  the  following  16  fundamental  sums  that  appear  above: 

T,  1  ?S  i 

stf.  fra1,  ?« , 

fc,£ ,  f  m , f  HT*. |S  s ,  ? IW«  > 

£  gp ,  2  iSTt'i;  5 ,  ILiif1,  2  IS  f,  f  li?  Ki;  S  5 ,  5  »l‘  (A-24) 
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It  should  be  observed  that  the  fundamental  statistical  quantities  required  of  random  variables  {zk} 
number  only  5,  namely 

(A-25) 

However,  they  combine  in  various  ways  to  yield  the  16  sums  defined  in  (A-24),  which  in  turn  are 
encountered  in  (A-19)-(A-23). 

In  terms  of  the  sums  defined  in  (A-24),  the  sum  of  (A-19)-(A-23)  can  be  written  as 

s  -  IfM*  •St**'-* 

V -  ♦  fr  l tb  “ 1 2  TJi  -  \  +  *  Xm  +  «  *  fcl  -  6  ^ 

++r4(v^-t»+2Th)^} 

+  2trM-Tj  +  |VT„r 

+  +(T„-Tj|T,r  +  2  HsfCV^TT'}  +  W.  (A-26) 

This  is  the  general  relation  for  the  average  of  the  magnitude-fourth  power  of  the  sum  (A- 1 8)  of  in¬ 
dependent  complex  random  variables. 

Special  cases  of  (A-26)  are  afforded  by 

$  *  \  *  W  for  N*  i  ,  (A-27) 

which  is  obviously  correct,  and  by 

5  »  V  2  T4J-  +  2T,’  +  IT,,)*  fw  %  .  o  4r  all  k,  (A-28) 

which  is  identical  to  (19)  in  reference  2  for  c  =  0  there.  That  is,  (A-28)  applies  to  zero-mean  random 
variables  {zk}.  If  we  let  complex  random  variable 

C  ■  ,  (A-29) 

then  (A-26)  is  obviously  an  expression  for  |  C|4.  At  the  same  time,  (9)  and  (A-24)  yield 

(A-30) 

Therefore  (5),  (A-26),  and  (A-30)  yield 

Vo  [ I  cfl  *  T4)  -  +  R.  fo]  -  2  Tw  -  T„  +  »  T.  +  t  R. ft)  -  (,  %, 

-1 

+  fe-T„j  4  |T„-  Tj’ 

4  2(t>.~T„)|t,)j  +  2ReI(T»-TH)Tf  J  (A-31) 

A-5 
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This  is  the  general  relation  for  the  variance  of  the  magnitude-squared  value  of  the  sum  (A-29)  of 
independent  complex  random  variables. 

Special  cases  of  (A-3 1 )  are  afforded  by 

Vo<£|C|l]  *  T4(-Tj  «  4r  M  =  l,  (A-32) 

which  is  obviously  correct,  and  by 

fer  ?*0  4  oil  k,  (A-33) 

which  is  identical  to  (22)  in  reference  2  for  c  =  0  there.  That  is,  (A-33)  applies  to  zero-mean  random 
variables  {zk}. 

A  program  that  evaluates  the  general  expression  for  the  variance  in  (A-31)  is  given  in  table  A-l .  The 
fundamental  input  statistics,  (A-25),  must  be  input  in  lines  90-190.  This  program  has  been  checked 
throughly  for 

1 .  N  =  1 ,  arbitrary  moments 

2.  N  =  2,  arbitrary  moments 

3.  N  =  3,  arbitrary  moments 

4.  arbitrary  N,  complex  Gaussian  random  variables  with  correlated  real  and  imaginary  parts,  and 
nonzero  means.  (See  appendix  B  for  this  derivation.) 

In  terms  of  the  fundamental  sums  defined  in  (A-24),  the  quantities  encountered  in  (7)-(ll)  are 
expressible  as 


Avfc}-  T, 

Me1}* 

Mic|‘)  =  t„-t„  +\t,)* 

Vor  {  C  j  »  l,  -  Tlt 

C/}  =  T„-  T„,  (A-34) 


A-6 
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Tabic  A*l.  Program  for  Variance  (A-31) 


ie  I  MOMENTS  OF  C  ■  SUM  OF  N  INDEPENDENT  COMPLEX  RANDOM  VMRIftBl.ES. 
20  I  LINES  50-190  GENERATE  INPUT  STATISTICS. 


30  OPTION  BASE  1 

40  DIM  Ztr<9>|21t(9>l Z2«(9 > , Z2r<9>,Z2>t9^, 
50  N>8 

S0  REDIM  Zlr<N>|Z|l<N>, Z2»( N) , 2Zr <H> , Z2 1  <M 
70  I  TO  INPUT  THE  FUNDAMENTAL  STATISTICAL 
00  I  OUANTITIES,  REPLACE  LINES  90-190. 

90  RANDOMIZE  SQRC.S) 

100  FOR  K<*1  TO  N 

110  ZIptK'-RND 

120  Zli<K)«RND 

130  22ft<K>*RND 

140  Z2r(IO-RND 

190  Z2i <K)*RND 

1E0  Z3r  < k  >*RND 

170  23t <K)*RND 

100  Z4<K>«RND 

190  NEXT  X 

200  FOR  K» 1  TO  N 

210  Ql*21r <K)A2 

220  02*21 1  <K>-'2 

230  Op*Ql+Q2 

240  Qft*Ql-Q2 

290  Qc*2*Zlr<K)«Zn  <K> 

200  Tlr«TlP*Zlp(K> 

270  T1 1  «T1 i *21 1  <K ' 

200  T21-T21*Z2»<tO 

290  T22p«T22p*Z2p<K> 

300  T221 “T22<  *Z21 <K> 

310  T23>T23*0p 

320  T24p»T24p*0» 

330  T241 »T24i +0c 

340  T31r»T31r*Z3r<K) 

390  T31»-T31l*Z3i<K> 

300  T32r*T32p+Z2«<K)»ZIrfK> 

370  T32t  *T32t ♦Z2ft\tO *21 1  <K > 

300  T33p»T33p*Z2p<K)*ZIp'K>*Z2i  <K  >»21  i  «.K> 

340  T33t»T33l*22i  <IO«Zlr  <K>-22p>  k  >*2 1  i  (K  > 

400  T34p*T34p+Qp*Z1p<IO 
410  T341 «T341 +Qp*Zl I  <IO 

420  T41*T4t*Z4<iO 

430  T42r»T42i*Z3r<K)«Zlr<K»*Z3t (F>*Z1 i (k> 

440  T43*T43+Z2»<K>-2 

490  T44-T44tZ2rfK)-'2+Z2i  <K>-2 

400  T49«T45+Z2»CK)*Qp 

470  T40r»T4Sr*Z2r<K>*Q»,*Z2i  <K>*Qc 

400  T47*T47+Qp*2 

490  NEXT  K 

500  Rl-TtP-2 

910  R2»T1W 

520  R3-T21-T23 

530  N4*T22p-T24r 

940  R5»T22I -T24 i 

550  Rp»Rl*R2 

500  R«*Rt-R2 

570  Rc«2*Tlr»Tlt 

500  V«T41-4*T42r-2»T43-T44*8»T45+4*T46r-6*T47 

590  V*V*4*<T31p-2*T32p-T33p*2«T34p>»Tlr 

000  V*V*4*<T31t-2»T32i-T33t*2*T34i >*T 1 1 

010  V*V*R3''2*R4A2*R9'2 

020  V-V*2*<R3*Rp4R4»P»+R5»Rc > 

030  PRINT  *N  -"JN 

040  PRINT  " Av<C  >  •  “ ;  Tlr;  “♦  i  <  “ ;  T1 1 J  “  > 

050  PRINT  "Av<CA2>  «  " JR4*R»;  “♦  t  <  “ ;R5*Rc  {  "  ‘ 

000  PRINT  "Av<C«4g-2)  »“IR3+Rp 

070  PRINT  "Std  D»u(C»»gA2>  ■"{SORiVj 

000  PRINT  "V*p<C»»gA2>  «"IV 

090  PRINT  "Av<C»4g'-1>  »"i <R3*Rp>- 2*V 

700  END 


N  •  8 

Av'C  «  3..  407426 1  3063  *U  4.64537563274  • 

Am*.  C'£  >  *  -4.57145748691  ♦(<  32.480628  3278 
AvCm  !<  *  33.4546032912 
Std  D*v<C«»(|'2>  *  6.7210029004 
V4rvCft*gA2>  *  45.1718799872 
AvfC»40'>4 >  »  1164.38236136 


3r<9..,23H9».Z4<9» 

,Z}p<H>,Z3t<N>,Z4<N» 
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APPENDIX  B 

COMPLEX  GAUSSIAN  RANDOM  VARIABLES 

For  complex  Gaussian  random  variables  {zk},  the  means,  variances,  and  covariances  add;  thus  all 
we  need  to  do  is  evaluate  these  quantities  for  a  typical  random  variable  z.  We  let  the  means  be  noted 
by 

+  »  «+  it.  (B-l) 

If  we  denote  the  zero-mean  components  of  the  real  and  imaginary  parts  as 


f  >a-k. 


then  the  variances  and  covariance  are 


V*  *  ,  f  »  (B-3) 

To  check  the  general  results  in  appendix  A,  we  need  the  five  statistical  quantities  listed  in  (A-25). 
They  are 

Avl*]  -  S  =  oM-it,  (B-4) 

.  a74  of  4  b*  +  .  (B-5) 

AvM  *  **  «  (x  +  ij)1  =  P-J  +.2*5  -  &+*)*- fry)* 

+  i2(a+*^(by)  =  a‘  +  (r*- (1?+^)  +  i2(oi  (B-6) 

fciV  *  ()l‘+31)(x+i^j  *  X5  4xjjr'  4-  iij3  (B-7) 

The  components  required  here  are  given  by 


X3  ■  fa  4-  *  Q*  4- 

-  Q  b1  4-  a  0-jj*  4-  2  y  <r«*j 

jjT-  b34  3b r* 

"ipy  s  cfb  4  br*  4  2a^>rnr3 . 


Av^W1*]  ■!*?>  *  £^4  ab’  4  a  (3t,J 4  2 *;  ** 

+  i£l?  4a'b  4  b  (r»4  3^)4  2a^  *v«ji] 


Hence 
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Lastly  we  need 


(*1*  »  (x* +•)’)*  *  ^  4-  y4  4-  2*Y , 

for  which  the  components  are 

^-£777-  a4  +  tirt+Zc* 

b4  +  t?*j+3$ 

»  {ft,4  2av+«‘Xl',+  2y  Y’) 

*  rftf  +  ciV*  +  +«;**,*  +'i^T*Cy. 


Then  there  follows 

Var[w*}  =  |?14  - 1*1*  -  2  (<r44<r4  +  yVfsj)  4-  *(aV* 4 ^ 4-y a«i  K) 

and 


fef  -  V{Wl]  +Vor?M’}. 


(B- 1 0) 


(B-l  1) 


(B-12) 


(B-l  3) 


A  program  that  incorporates  the  above  relations  is  presented  below  in  table  B-l .  The  Gaussian  rules 
occur  in  lines  280-320. 


aiwwywff  -y  whmr  11  1 eiiiv^  ymywaigjTWIiragp^BPIH^ 
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Table  B>1.  Check  Program  for  Complex  Gaussian  Random  Variable! 


10  I  MOMENTS  OF  C  •  SUM  OF  N  INBEPENBENT  COMPLEX  PANl-OM  VARIABLES. 
20  i  GAUSSIAN  TEST  EXAMPLE. 

30  OPTION  BASE  1 

40  Bill  2lr<9),Zli  ^S  > ,  Z2«t<9> ,  22r  <9> ,  22i  *' 9  >  .  Z3.'«9>  ,231  <9>,24(9> 

30  N*8 

60  REBIM  ZIr<N>,2Ii<N>,22»<N),22r(N>,22i  <:u  >  tZ%r  <  N.\  Z3i 'N> ,  24<N> 

70  RANDOMIZE  SQR(.6> 

80  FOR  Ml  TO  N 

90  A-2+RNB-1 

100  B«2»RND* 1 

110  Si gx»RNB 

120  Stgy«RND 

130  Rho«2*RNB~ 1 

140  A2«A-'2 

ISO  B2«B2 

160  Vx«Sigx'2 

170  Vy«SigyA2 

180  P«Rhc>*Sig><#Stgy 

190  Zlr <K>»A 

208  21 l <K>»B 

218  Z2*i<|0*>fl2+Vx+82+Vy 

228  22r(:tO*>A2+Vy".B2+Vy:< 

230  Z2i  <lO-2«<A*F  +  P> 

240  23r <K)*A»<A2+B2 >+H»< 3»Vx+Vy >+2*B*P 

250  23 i  <K>»B  +  <A2  +  B2>+8(K  Vx+3»Vy  :■ +2*  A+F' 

260  V*r»2+C  VxA2+Vy ''2+2»P,'2>+4*'.  A2*Vx+B2*Vy*C  +  P«-8*P  ‘ 


278 

24  tk i*Z2ni<K>'2+V*r 

280 

G»G+fl  ! 

GAUSSIAN 

RULES 

290 

H«H+B  ’ 

GAUSSIAN 

RULES 

300 

Vd-Vd+V;  i 

GAUSSIAN 

RULES 

310 

V»»V*+Vy  • 

GAUSSIAN 

RULES 

320 

Vc-Vc+P  ! 

GAUSSIAN 

RULES 

330 

NEXT  K 

340 

FOR  Ml  TO  rl 

350 

Ql»Zlr<tOA2 

360 

02-21 1  <K>'2 

370 

0p-Q  1  +02 

380 

Qm»01-02 

390 

0c«2*Zlr  <10*21 1  <►> 

400 

Tlr<iT  t  r+ZIr <k  > 

410 

Til  »T  l  ‘  +21 1  <IO 

420 

T21  »T21  ♦?2ii»<k> 

430 

T22r»T22r+Z2,''K> 

440 

T22i»T22i+Z2t  ' k ' 

430 

T23*T23+0p 

460 

T24r«T24r +Um 

470 

T24i-T24i+Qc 

480 

T31r-*T31r+23r<k> 

490 

T31 i«T31 i+23i <k> 

500 

T32r  »T32r +Z2(ii  <k>*21r(K 

) 

510 

T32i«T32i  +Z2m<M+Zl i 'K> 

320  T33r«T33r+Z2r'<K.+Zlr<:iO+Z2i  <F.)*2I  i  <K> 

330  T33i*T33i +221  00#21r<k  <-Z2r < K >*Z1 1  >K> 

340  T34r«T34r+Qp*Zlr<K> 

530  T34i *T34i +0p»Zl i <K> 

360  T41-T41+Z4<K> 

570  T42r«T42r+Z3rtlO*Zlr<K>+Z3i <K>*Z1 i <K> 

380  T43»T43+Z2(u<K  >A2 

590  T44«T44+22r<k>'2+Z2<  >  tO-2 

600  T45«T45»22i»<:K>*Qp 

610  T46r»T46r+Z2r<K)»0«i+22(  UO*Oc 

620  T47»T47+Gp  2 

630  NEXT  6 

640  Rl*Tlr  2 

650  F!2*T  1  i  "2 


pn  iwu'wi  iimump— 
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Table  B-l.  (Coat'd)  Clack  Program  for  Complex  Gaomlaa  Raadom  Variable! 


670  R4«T22r-T24r 

660  R3»T22l-T:J4i 

see  ftp*Ri+R2 
roe  r»»»ri-ri 

710.  Rc-2«Tlr*Tli 

720  V»T41-4»T42r-2*T43-T44*8#T45<-4*T46r-6*T47 

730  V«V+4*<T31r-2*T32r-T33»'+2*T34i'>»Tlr 

740  v»V44*<T3H-2»T32i-T33i+2»T34i>»T11 
7S0  V«V*R3*2*R4- 2+R3 '2 

760  V*V+2*<,R3*Rp*R4*Rfn+R3»fic ) 

770  PRINT  "N  »“{N 

780  PRINT  "flv<0  •  "  J  T  l r ;  11  ♦ )  i.  ”  ;  T 1 1 J  "  ■<  ” 

790  PRINT  “flv(C-'2:'  «  "  J  R4+Rmi  "  +  i  <  "  J  R5*Rc  {  "  ■ 

800  PRINT  "A<.'<Cfn4gA2>  ■  ";R3VRp 
810  PRINT  "Sid.  D*w'Cm*g  2>  »";S0R<:V> 

820  PRINT  "V*i',C«.«9  2)  «"}V 

830  PRINT  v4)  «“} <R3+Rp>'2+V 

840  PPINT 

830  PRINT  "USING  GAUSSIAN  RULES:* 

860  PRINT  Gl“*1<";H;">" 

870  PRINT  G  '£*Vd-<H-'2+V*:>;  "  +  t  < '  ;2*<G»H+Vc  '  V 

880  PRINT  GA2+Vd*HA24 V4 

890  V*rc»2*'  Vd'2+V«~2+2*Ve  2^*4*<G-'£*Vd  +  H  £*V«+2*G*H»Vc  > 

900  PRINT  ?0R<V»rc> 

910  PRINT  V»rc 

920  PRINT  <G'  2+Vd+H  2+V4>'-2+V*cc 
930  END 


M  «  8 

Av t C )  «  -.81709993443  +  i<  .81893514539  > 
P«‘<C  2>,»  -.937717628742  ♦  t  < -.  76 1  38494979  \> 
flv"'Cni*g  2)  »  6.  17375336444 
Std  D*v<Cm*g-2>  *  6.80940237934 
V*r7C*i»g-'2>  *  36.  1 129137351 
fl'.'(Cm»g''4 >  *  74.2523433733 

USING  GAUSSIAN  RULES: 

-.81709993443  t-U  .81395514559  > 

-. 98771762373  ♦ t 761384949773  • 

6. 1 75 75336445 
6.00946227937 
36. 1129157554 
74.2528453739 
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APPENDIX  C 

EFFECTS  OF  POSITION  AND  DELAY  PERTURBATIONS 
ON  LINEAR  ARRAY  RESPONSE 

Consider  a  plane-wave  arrival,  at  frequency  f„  from  angle  propagating  with  speed  c,  as  shown 
in  figure  C-l.  If  the  propagation  time  delay  to  reach  the  origin  O  is  defined  as  aero,  then  the 
propagation  delay  to  reach  position  x  on  the  linear  array  is 

TM  (c-i) 


Figure  C-l.  Linear  Array,  not  Necessarily  Equispaced 


The  processing  of  interest  is  depicted  in  figure  C-2.  The  position  of  the  Jc-th  element  in  the  linear 
array  is  xk;  the  delay  employed  in  the  k-th  branch  of  the  beamformer  is  Dk  (0<),  which  depends  on  the 
desired  look  direction  0ti  and  the  corresponding  element  weight  is  wk.  The  voltage  transfer  function 
of  the  propagation  delay  and  beamformer,  as  applied  to  the  arriving  plane-wave,  is  then 

V«  (C-2) 


BEAMFORMER 


Figure  C-2.  Processing  of  Plane-Wave  Arrival 
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Now  suppose  that  the  desired  position  of  the  k-th  element  in  the  linear  array  is  dk,  but  that  the 
actual  position  is 


X*  -  4  ♦  a* ,  (C-3) 

where  Ak  is  a  random  (distance)  variable.  Notice  that  we  are  restricting  consideration  to  longitudinal 
element  movement  only  along  the  linear  array;  this  restriction  is  removed  later. 

Also  suppose  that  the  desired  k-th  element  time  delay  to  steer  the  array  in  direction  0,  is  (ignoring  a 
common  bulk  delay  for  all  elements) 


-  sin  ^  ,  (C-4) 

but  that  the  actually  employed  delay  is 

JKW  ■  -  4->'»  S  *  5.  ,  <c-5> 

where  dk  is  a  random  (time)  variable.  (Notice  that  this  choice  of  element  delay  takes  no  account  of 
distance  perturbation  Ak;  a  case  where  this  distance  perturbation  is  taken  into  account  is  considered 
later.) 

The  transfer  function  of  the  linear  array  is  then* 

V,  *  Sh  «r[-  i  2H;  w  ft  f  i  2H  $  -  i  H  $ ,] 

‘  [l  2rrf.  -  i  J'i»  ft,  +  S,)]  (C-6) 

The  effect  of  the  position  and  delay  perturbations  leads  to  the  k-th  random  phase  shift  term  exp(i<)k], 
where 

*  -  2irf,  JW  ft,  +  -  ”  Ar  (-^  SM  ft  +  -^).  (C-7) 

Here  A,  is  the  wavelength  of  the  plane-wave  arrival,  and  Ta  is  the  period  of  the  arrival.  The  fraction  of 
a  wavelength  movement,  Ak/Aa,  and  the  fraction  of  a  period  delay,  dk/Ta,  are  obviously  important 
parameters. 

Now  if  the  position  and  delay  perturbations  have  the  properties  t 

bl\)  =  o,  DwfAi]  -  oj  1  f  „ 

TW  a);  k, 

Av  ^4  =  o ,  5U  Dc*  m  *  (C-8) 

✓ 

and  Ak  is  uncorrelated  with  dk,  then  we  find 


*  Weight  perturbations  and  element  failures  are  easily  incorporated  by  replacing  wk  by 
wk(l  +  rk)gk,  where  rk  is  a  zero-mean  random  variable  and  gk  is  a  (0,  1)  random  variable;  these 
features  are  not  included  in  this  appendix, 
t  These  could  be  generalized  to  allow  dependence  on  k. 
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Avft)  •  0,  Sto  bwfo  *  2ir  [fj-j  4^.J]  *  ^  ,  (C-9) 

which  is  a  function  of  arrival  angle  6t.  Thus,  if  we  plot  average  properties  of  the  beamformer,  by 
varying  the  look  angle  9t  while  holding  arrival  angle  0a  fixed,  we  can  treat  o^asa  constant.  However, 
if  we  do  the  reverse,  we  must  take  into  account  a  varying  with  arrival  angle.  Of  course,  if  od  ■  0  (no 
position  perturbations),  then  o+  -  2n  o,/T,  is  a  constant  independent  of  any  angles.  We  might  also 
note  that  for  broadside  arrivals,  9a  -  0,  o+  is  a  minimum  and,  in  fact,  positional  perturbations 
(od  *  0)  do  not  affect  the  transfer  function,  (C-6). 

Instead  of  choosing  element  delay  Dk(6|)  without  taking  account  of  position  perturbation  Ak, 
assume  that  we  have  this  information  (perhaps  via  measurements  after  construction  of  each  array) 
and  decide  to  use  it  by  employing,  instead  of  (C-S),  the  time  delay 

&(<*)«- ^**6, +  (C-10) 

where  dk  is  again  a  random  (time)  variable.  (Of  course  { Ak}  could  still  be  random  variables  from  one 
array  construction  to  another.)  Then  the  transfer  function  of  the  linear  array  is 

^  *  SF*  '2*4.  4  i^irt $ -  i  2*4.  i,"] 

»  txp[i2ir^%sVV ftr3Hl4tyl  (C-ll) 

The  effect  of  the  position  and  delay  perturbations  then  leads  to  the  k-th  random  phase  shift  term 
exp[i+k],  where  now 

"  2irf,  +  £)  =  *  (l^  [/>"  ft,'  .  (C- 1 2) 

Keeping  the  same  statistical  properties  as  assumed  above,  we  now  find 

Sk  Dwffc}*  ^ +  (vJ]  *  ,  (C-i3) 

which  is  now  a  function  of  both  arrival  angle  0a  and  look  direction  9t.  If  od  «  0  (no  position  per¬ 
turbations),  then  0|  -  2n  ot/Ta  is  independent  of  any  angles.  For  od  ¥•  0,  the  minimum  of  o+  is  realized 
when  6,  -  0a . 

For  positional  perturbations  into  the  x,y  plane  (rather  than  just  along  the  x-axis),  we  have  instead 
of  (C-l),  propagation  delay 


~~i  -  \  >!  S**Qa  C»S  fti 

T(*,  U,  ft,)  =  - ^ -  (C-14) 

Then  the  voltage  transfer  function  of  the  linear  array  becomes 

V  =  -Sf*.  [-  i  M  t(\  * ,  ft,)  -  I  2trf,  X\  (e,)]  .  (C-l  5) 

Now,  instead  of  (C-3),  let  element  position 

\ *  N  ,  (C-16) 


C-3 
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where  Ak  end  *,  are  random  (di.tance)  variable,.  Then,  If  we  still  use  the  ..me  delay  a,  In  <C-», 

4l*'"  ^  ■  (c-i7) 


we  find  the  perturbed  transfer  function  to  be 

%  „  ^  7^ *■(*'»  4' 3"  *)]  **f 040  , 


where  now 


^  -  M.  _  4  **  ^  ^ 


(C-18) 


(C-19) 


Analogous  to  (C-8).  if  we  now  assume 

x  r  „  ~r  «-■  T*1.  r*  T*.^  aX  * /°»r» »  (C-20) 

5^,0, 0,V0,  ^y  th  -  V  H  f 

where  nonrero  „  corresponds  to  correlated  x,y  movement  of  an  element,  we  find 

AvHI)  -o, 

A 

.  (C-21) 

This  is  independent  of  «„  as  was  (C-9);  thus  the  comments  directly  under  ,C-9)  me  relevant  in  this 


(C-21) 
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APPENDIX  D 

VARIANCE  OF  POWER  RESPONSE 

LINEAR  ARRAY 

The  variance  of  the  power  response  was  given  in  (A-31).  We  now  specialize  it  to  the  equispaced 
linear  array.  Define 

C,  *  »  C**  “  » 

C3  >  Q*  (D-l) 

Then  the  various  quantities  {Tjj}  in  (A-24)  may,  with  the  aid  of  (18)  and  (40)  and  the  realness  of  y j »/2» 
be  expressed  as 

T,  *  S  l, 

C,*  K  ,  T„  a  cir  L,  ,  T.J  *  ClK,  T*-  c*  L. 

T„  *  Cj  Lj  ,  Tj,  r  C,  Tj,*  C,  C„  Lj ,  V  -  Cf 

t4(>  Qw4,  t*  •  V  CK,  l-c3;w4 

T„  -  c *  UK ,  \  -  c; c,r  w4 ,  r„ «  cf  w4  (D-2) 

The  only  terms  in  (D-2)  that  depend  on  u,  areT,,  T22,  T24,  and  {T,j}f,  by  virtue  of  depending  on  L,, 
L2,  Lj.  Now  (A-3 1 )  can  be  expressed  as 

VarilCl*}  =  V  4T41-2T4J-l^f  ?r4S4tTH-fcT„ 

♦  (r„  Tj  ♦  iTu - tJ *  z (t t„-  r„)  r; 

*  (V  cj  +  sc; c2.+-4c,VM-  icOk 

+  4  C,  (v,  -  2C,  Clm  ~  C,  C»r  +  2  C,J)  L,  L , 

+  (u-  c?JV  +  (u-  c*)‘  L* 

+  2  c“f[(U-  c;)W.f  (c„-  C,0lJ  L?  •  (D-3) 

in  terms  of  the  quantities  in  (D-l )  and  (40).  Finally  we  rearrange  (D-3)  to  read 

VarJjcrt  *  F,  L,l,  +  f,.  l,‘  ♦  F„mu:  ♦  F*  l,  L|  ,  (D-4) 


D-l 


mmw  ip  .1  wwiiMwaiiP  u  »i«w*  nwuvmpw 


Mpi««il!iuiii  \  w 
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where 

F4  »  C,  -  *C,  C»  -  2  C*  -  C,:  V  ?  Cfc,.  +  c,r -  6  c r 

F  •  VC,  {Cr2CSn»-  C,C*r+2C?) 

Fa«-  fc-  C0’  ,  Fu-  fcr-  C$ 

F*  •  *  crfc.-  C*)  ,  Fk  -  2  cr  (Q,  -  c;)  (D-5) 

Since  L,  and  L3  only  occur  in  the  combinations  Lt  L3  and  Lf,  we  note,  using  the  properties  of  L,,  LJt 
L3,  that  (D-4)  is  even  in  u  and  has  period  2n  in  u.  Therefore  we  need  only  to  evaluate  (D-4)  in  the  range 
(0,  n)  for  u  for  this  equispaced  linear  array.  The  program  for  evaluating  (42)  and  (D-4)  is  presented 
below. 

DEEP  SIDELOBE  REGION 

In  the  deep  sidelobe  region,  L,,  L2,  and  L}  of  (40)  are  approximately  zero;  then  the  only  {T^}  in  (D- 
2)  that  are  nonzero  are  T2I ,  T23,  and  { T4J } .  It  then  follows  from  (D-4)  that 

Varllcr]  *  W4  +•  Wi  in  si  JtloifS,  (D-6) 

where  these  quantities  are  given  in  (D-5)  and  (40).  The  constants  needed  in  (D-5)  are  given  by  (D-l). 
No  simpler  expression  for  (D-6)  appears  possible. 

Alternative  expressions  and  interpretations  for  F2a  and  F4  in  (D-6)  are  possible;  let 

-2*  *  p*  »  p.  *  (*+  ri)  0  O  >  (D-7) 

where  we  used  (16).  Also  let 

V  P«-  '  (D-8) 

Then  we  find 

Fa.  -  Kf  ,  5  *  Kf  -  2  Kf  -  .  (D-9) 

Comparing  these  results  with  (40)  and  (44)  in  reference  2,  we  see  that  and  F4  are  equal  to  the 
cumulants  A, ,  and  A22  of  random  variable  pk. 

PLANAR  ARRAY 

The  variance  of  the  power  response  is  given  in  (A-31);  we  now  specialize  it  to  the  planar  array.  With 
the  aid  of  (A-24),  (18),  (D-l),  and  (75)  and  the  realness  of  yM  y2,  there  follows 


D-2 


*m  i  mi/m  iphi  !^h  'P” 
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i;  -  C,  l*Im>  l%) 
r„-  C^W,*,  T„-C„  L*WiSfW 

rn  •  c,W,  t„  *  c,‘  uwi»»lTW 

[T,))  {C,,f,C,,,C,Clr,  C,*j L* Wl? W 

W’  •fc.c.S.C.C, cx.cfe,  c fl'tf'MC  .  (D.10> 

The  terms  that  depend  on  u  and  v  are  T,,  T2J,  T24,  and  {T3j}f,  by  virtue  of  depending  on  L,,  L2,  L3, 
defined  in  (75). 

Now  (A-31)  can  be  expressed  as  the  upper  half  of  (D-3),  which,  in  turn,  can  be  developed,  as  above, 
into 


Vorild'j  .  FXV+  f,  l7(m)l“wl“wl*(v) 

t  WWf’  +  F»Lf(u)LfW 

F»  L*l“)L?(>j Lf (u)C'M  ,  (D.,1) 

where  the  F-constants  are  given  by  (D-5).  This  is  the  final  result,  where  the  necessary  functions  are 
defined  in  (75). 

Programs  for  linear  array  results,  (42)  and  (D-4),  and  for  planar  array  results,  (76)  and  (D-ll),  are 
presented  below  in  tables  0-1  and  D-2,  respectively.  For  the  linear  array  case  in  table  D-l,  the  inputs 
required  of  the  user  are 

H  in  line  20,  the  number  of  elements  on  one-half  of  the  linear  array 
Q  in  line  40,  the  probability  of  element  failure 

Sigmar  (or)  in  line  50,  the  standard  deviation  of  the  relative  error  of  the  weights 
Sigmap  (Op)  in  line  60,  the  standard  deviation  of  the  phase  perturbation  in  radians 
{ wk}f*  in  lines  80-90,  the  H  weights  on  one  half  of  the  array. 

Alternative  weight  structures  are  available  in  lines  220-240,  if  desired. 

For  the  planar  array  in  table  D-2,  the  inputs  required  of  the  user  are 

Phi _ t  (ft)  in  line  20,  the  polar  look  angle  in  radians 

Theta _ I  (0f)  in  line  30,  the  azimuthal  look  angle  in  radians 

Slice  =  1  or  2  in  line  40,  depending  on  polar  or  azimuthal  slice 

Theta _ a  (0,)  in  line  60,  the  azimuthal  arrival  angle  or  Phi _ a  ($,)  in  line  80,  the  polar  arrival  angle 

Hx,  Hy  in  lines  90  and  1 10,  the  number  of  elements  in  the  halves  of  the  array  in  the  x  and  y  coor¬ 
dinates 

Dx _ lam,  Dy_iam  in  lines  130  and  140,  the  ratio  of  x  and  y  array  spacing  to  the  arrival  wavelength 

{wJ*)}Fx,  {w^>}py  in  lines  150-200,  the  weight  structures  in  the  x  and  y  coordinates 
Q  in  line  210,  the  probability  of  element  failure 

Sigmar  (or)  in  line  220,  the  standard  deviation  of  the  relative  error  of  the  weights 
Sigmap  (Op)  in  line  230,  the  standard  deviation  of  the  phase  perturbations  in  radians. 
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T*bk  D-l.  Program  for  Equlspaced  Linear  Amy 


to 
£0 
30 
40 
30 
<50 
70 
30 
30 
100 
no 
1 20 
130 
140 
ISO 
160 
170 
ISO 
190 
*00 
*10 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
4~0 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 


•  RANDOM  BEAMFORMER;  EQUI-SPACED  LINE  ARRAY 
H» 1 O  !  HUMBER  OF  ELEMENTS  ON  ONE  HALF  OF  ARRAY 

!  SYMMETRIC  WEIGHTS;  TOTAL  NUMBER  OF  ELEMENTS  «  3H 
Q-.OOl  (PROBABILITY  OF  ELEMENT  FAILURE 

Si  gm  *1  -.01  I  STANDARD  DEVIATION  OF  RELATIVE  ERROR  OF  WEIGHT 

Sigm»p«.01  I  STANDARD  DEVIATION  OF  PHASE  <  PAD  1 A  NS  > 

I  DOLPH-CHEBYSHfev  WEIGHTS  FOR  -30  DB  SIDEL0BES1 
DATA  2.405798938,2.333866168,2. 195114975, 1 . 999277697, 1 . 759768704 
DATA  1.492415162, 1.213997044, .9407573274, .6870419302, . 7833503571 
PRINTER  IS  0 

PRINT  "EQUI-SPACED  LINE  ARRAY  JOTAL  HUMBER  OF  ELEMENTS  ■" 

PRINT  "Probability  of  element  failure  «";0 

PRINT  "Standard  deviation  of  relative  irror  of  *>*  i  ght  “"(Sigmar 
PPIHT  "Standard  deviation  of  phase  <tn  radians)  ““;Sigmap 
DIM  Mu< 1 14) ,Nu< 1 :4>,Cammat 1 : 2) , Nt 1 : 50) 

DIM  Ideal  power  <. O:  256) ,  Auer agepower  (0:256) ,  Si  gmapower<0:  256) 

Num«256  I  NUMBER  OF  PLOTTED  POINTS  O  256 
REDIM  WUIH) 

READ  W<*> 

T»0 

FOR  Ml  TO  H 

!  M<K  >“. 34  +  , 46*C0S(P1*(K-.  5)  •  H)  !  HAMMING  WEIGHTS 
!  MOOa  t  »COS<  P I  5)rH)  !  HANNING  WEIGHTS 

l  W ( K > *  1  1  UNIFORM  WEIGHTS 

T  =  T  +  U<.  K ) 

NEXT  K 
W2»W4«0 
FOR  Ml  TO  H 
W  <  K  )  ■  W  <  K )  r  T 
F»H< K  )  '-2 
W2*W2+P 
W4«W4+PA2 
NEXT  K 
W2»W2/2 

PRINT  "Effective  number  of  elements  ;  1  M2 
PRINT  LINY  1 ) 

W4«W4.-3 
Vr*S i gmar  a2 
Vp»S  i  g<nap'-2 

Mu  <  1  >  •■’Mu  <  2  )  “Mu  <  3  >  "Mu <  4  )  “  1  -0 

Hu i l'a 1 

Nu ( 2  >  *  1 >Vr 

Nu<  3 ’ « 1 +34Vr 

Nut 4 ) ■ 1 +6* Vr +  3*Vr  2 

Gammat 1 )*EXP<-. 5*Vp) 

Gamma<  2 ) “EXP( -2*Vp ) 

Cl ■Mu< 1 >*Nu( 1 )*G*mma( 1 ) 

C2m“Mu<2>*Hu> 2) 

C2r“M'i(  2  )*Nu<2)  *  Gam  mat  2 ) 

C3«Mu<3)*Nu(3>»Gamma( 1 ) 

C4*Hu<  4 ) *Nu<4 ) 

C 1 2“C 1 A2 
Co“2*C2»>tC2r 

F4»C4~4»Cl*C3-2*C2m''2-C2rA2+4*C12*Co-6»|.  1*2 
F3“4«Ct*<C3“Cl*Co  +  2*C12*Cl  1 
Aa®C2m-C 1 2 
Ab-C2r-C12 


F2a“A*-2 
F2b‘ Ab  * 
Fla“2*C12*Aa 
Flb“2*C12*Ab 
Com«F'4*W4 tF2a*rW2  2 
F 1 aw2*F 1 a*W2 
Var_vol t  *Aa#W2 
P»PT  burn 
FOR  I“0  TO  Hum 
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T«bk  D-l.  (Coal’d)  Program  for  EqvtspacMl  Llorar  Array 


6?o  u-i*p 

680  L 1  -L2*L3,*0 
600  FOR  K- 1  TO  H 
700  8»W<K> 

710  T-$*CQS<<.K-.3>*U> 

720  S-S->2 

730  Ll-Ll+T 

740  L2-L2+2*T-'2-$ 

750  L3-L3+S*T 

760  NEXT  K 

770  L2-L2--2 

780  L3«L3'4 

790  Id*«1 pow*r < I >-L1a2 

800  Rwtr»o*pow#r< 1 >»C12*Id*»1 pow«r< I )+V*r_/oH 

810  v*r  powtr*CoM+F2b*L2A2+F3*Ll*L3T<Fl*w2«'Flb»L2>*l 1-2 

820  Sign*pow«r<I)*SQR<V*r  pou*r> 

830  NEXT  I 

840  PLOTTER  IS  "GRAPHICS" 

830  GRAPHICS 

860  SCALE  0, Nun, -7,0 

670  GRID  Num/10,1 

880  PENUP 

890  FOR  1*0  TO  Nun 

90O  A-Idt«1pow*r< I > 

918  IF  AME-7  THEN  940 
920  PENUP 
930  GOTO  950 
940  PLOT  I , LGT <A> 

930  NEXT  I 

960  PENUP 

970  FOR  1-0  TO  Nun 

980  PLCT  I ,  LGT<Autr'*g*pow«t'<  l  /'> 

990  NEXT  I 

1000  PENUP 

1010  FOR  I«*l  TO  Nun 

1020  PLOT  I ,  LGT <Au*r«g*pow«r < I >+Si gm*pow»r < I  1 ) 

1030  NEXT  I 

1040  PENUP 

1030  FOR  1*0  TO  Nun 

I860  PLOT  I  ,LGT<Ruer*g*pow*r  <  I  >»2*$i  gn»pow«r  <:  I  >  > 

1070  NEXT  I 

1888  PENUP 

1090  PAUSE 

1100  DUMP  GRAPHICS 

1110  PRINT  L!N<8> 

1120  PRINTER  IS  16 
1130  END 


*  "T'T*1  »l»'.  ''f?*y).'r*n^^n!^m«ij''f  *«■ 
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Table  D-2.  Program  for  Equbpactd  Planar  Array 


10  1  PflNEiOM  BEAMFOPMER;  EUUI-SPACED  PLANAR  ARRAY 

£0  INPUT  "Polar  Look  Angle  <  n  radian**  ■  , Fni _1 

30  INPUT  "Azimuthal  Look  Angle  (in  radian;  >  *  ,",Tln.i  a  1 

40  INPUT  "  For  Polar  Site#  in  Arrival  Angle,  enter  7 

For  Azimuthal  Slice  in  Arrival  Angle,  enter  2“, Slice 

30  ON  Slice  GOTO  €0,30 

60  INPUT  "Fixed  Azimuthal  Arrival  Angle  fin  radian;*  »  ?",Theta  a 
70  GOTO  90  ' 

90  INPUT  "Fixed  Polar  Arrival  Angle  (in  radian**  *  ?",Pht  a 
90  Hx« 10  !  NUMBER  OF  ELEMENTS  ON  ONE  HALF  OF  ARRAY,  IN  X -DIRECT  ION 

100  !  SYMMETRIC  WEIGHTS;  TOTAL  NUMBER  OF  ELEMENTS  IN  X  -  2  Hx 

110  Hy«9  !  NUMBER  OF  ELEMENTS  ON  ONE  HALF  OF  ARRAY ,  IN  Y-DIRECTION 

1£0  !  SYMMETRIC  WEIGHTS;  TOTAL  NUMBER  OF  ELEMENTS  IN  Y  -  2  Hy 

130  Dx_  1  *M».  3  I  X-SPAC INGxWAVELENGTH 
140  Di*Jam«,5  !  Y-SPAC I NG/WAVELENGTH 

130  1  10  DOLPH-CHEBYSHEV  WEIGHTS  FOR  -30  DB  .SIDEi  OBES: 

ISO  DATA  2.405798938,2. 333866160, 2. 193114975, 1 . 939277697, 1 . 739768784 
170  DATA  1.492415162, 1.213997044, .9407373274, .6870419302, . 7833303571 
180  !  9  DOLPH-CHEBYSHEV  WEIGHTS  FOR  -30  DB  SIDELOBES! 

190  DATA  2.681 134180,2. 381677308, 2,391494643.2. 127093009, 1.810713745 

£00  DATA  1.467761548, 1.  123942032, .802315104,  .825036m87 

£10  Q«.01  !  PROBABILITY  OF  ELEMENT  FAILURE 

220  S i gmar * . 02  ISTANDARD  DEVIATION  OF  RELATIVE  ERROR  OF  WEIGHT 

£30  Si  ginap*.  1  ISTANDARD  DEVIATION  OF  PHASE  (.RADIANS* 

240  !  PRINTER  IS  0 

£50  PRINT  "Equi-»p*c»d  planar  array  Total  number  of  element*  ■  " ; 2*Hx; 2*Hy 

£60  PRINT  "X-spac  ing/wave length  »" ;  Dx_1  am,  "Y-sp*c  i  ngruavel  ergth  »"JDy_iam 
270  PRINT  "Probability  of  element  failure  »";G 

£80  PRINT  "Standard  deviation  of  relative  error  of  me i ght  « " ;  S  i  gmar 

£90  PRINT  "Standard  deviation  of  phaee  (in  radian;*  «";Sigmap 

300  DIM  Mu <1 1 4> , Nu( 1 ! 4 > , G amm  a( 1  i 2  > , Hx ( 1 : 50  > , Wv ( 1 : 50 • , U( 8: 256 > , V ( 0 ! 256 * 

310  DIM  Ideal  power (  0! 256* , Aver  agepewer f O: 256  * , Si gmapower ( 0 : 256* 

320  Num-256  i  NUMBER  OF  PLOTTED  POINTS  <=  256 

330  REDIM  Wx< 1 ! Hx* , Wv( 1 : Hy > 

340  READ  Wx  <  *  * , Wy  ( # ) 

350  T»1  ■/SUM(Wx  > 

360  W2.x  =  W4xa0 
370  FOR  K«1  TO  Hx 
380  Wx(K)“Wx<  K>*T 
390  P » W x •" K  * ' •  2 
400  W2x»W2x+P 
410  W4x»W4x+P  ''  2 
420  NEXT  K 
430  W2x»W2x/2 

440  W4x«W4x."8 

450  T*l/SUM(Wv'< 

460  W£y«W4y»0 

,70  FOR  (>1  TO  Hy 
480  Wv(K**Wy(K>*T 
490  P*Wv(K>A2 
300  W£y»W2y+P 
310  W4y«W4y+P-2 

320  NEXT  X 
330  W2y*M2y ' 2 
540  W4y*U4yr8 

350  PRINT  "Effective  number*  of  element*  in  x,v  - ;  l  W2;  ;  1  'W2y 
360  PRINT  LIMCl* 

370  PRINT  "Polar  Look  Angle  »";Phi  1 
580  PRINT  "Azimuthal  Look  Angle  «"7Th«ta  1 
590  ON  Slice  GOTO  600,630 

600  PRINT  "Fixed  Azimuthal  Arrival  Angle  =; Theta* 

610  PRINT  "Polar  31  ice  in  Arrival  Angle  from  '?  *.  PI  *" 

620  GOTO  650 

630  PRINT  "Fixed  Polar  Arrival  Angle  «";Phi_a 

640  PRINT  'Azimuthal  Slice  in  Arrival  Angle  from  u  to  2*PI " 

630  Vr “Si gmar ' 2 
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Table  D*2.  (Coat’d)  Program  for  Etgulspaced  Planar  Array 


660  Vp*Stgm*pA2 

670  MuC  1  >*MuC2>«MuC3>*MuC4>*l-Q 

600  NuCl>*l 

690  Nu<2i«t+\V 

700  NuC3 >*la3*Vr 

710  NuC4>*l+6*Vr+3*Vr "2 

720  Gamm*Cl>*EXPC-.5*Vp:> 

730  GammaC2>«£XPC-2*Vp> 

740  C 1  *Mu(l  'HrNu(  1  >*C*mm*C  1 ) 

750  C2m*Mu(2)*Nu<2) 

760  C2r«MuC2>*HuC2>4Camm*C2> 

770  C3*MuC  3)*Hu(3>*Cmm(  1  ) 

7S0  C4*Mu<4>*NuC4> 

790  C 1 2*C 1 A2 
800  Co*2*C2m+C2r 

810  F4*C4-4*C1  *C3-2»C2mA£-C2rA2*4*C  12*Co-6*'i  1<2  2 

820  F3*4*Cl*CC3-Cl*Co+£*C12*Cl > 

830  Aa-C2m~C12 
840  Rb«C2r-C 1 2 
858  F2*»AaA2 
860  F2b»AbA2 
870  F 1 a»2*C 1 2*8* 

880  F 1 bc2*C 1 2*Ab 

890  Com*F4*W4x*U4 y +  F 2** C  W2;<*W2y  >  A2 
900  F 1 aw2*F 1 a*W2x*W2y 
910  Var  wol  t*Aa#H2v,»W2y 
920  SI »SlNCPhi_l ) 

930  Ax*2*PI*0x~! am 

940  Av«2»PI  *Dy~l  am 
950  Ul«Ax*Sl*COSCTh*ta  1 > 

960  V1*Ay*S1*SINCThe<4~l > 

970  ON  Slice  GOTO  Pol ar, Azi nut hal 
980  Polar:  Ua«Ax*C03cTheta_a> 

990  Va*Ay*$IN(Th*ta  a> 

100O  P=.5*PI-Hum 
1010  FOP  I«0  TO  Hum 
1020  Phi  a»P*I 
1030  3a«SIHcPhi  a> 

1040  Ik  I  )*U1 -Ua*Sa 

1050  VC  I >*V1 -Va*Sa 

1060  NEXT  I 

1070  GOTO  Common 

1O80  Azimuthal:  Sa*SINCPht_a) 

1098  Ua«A**Sa 

1100  Va»Ay*Sa 

1110  P*2*P I  Hum 

ll20  FOR  I«0  TO  Hum 

1133  Theta  a»P*I 

1140  Ucl )«Ul-Ua*C0S' Theta  a» 

1150  Vc I >»V1 - V**S1H' Theta"*) 

1160  NEXT  1 

1170  CommonlFUR  1*8  TO  Hum 
1180  i. lu*L2u*L3u*0 
1190  Ui  *U<  I  ) 

1200  FOR  M*l  TO  H,' 

1210  S«W;;CK> 

1220  T*S*C0SccK-.5-*Ui ) 

1230  S*S  •  2 
1740  Llu-Llu+T 
1258  L2u«L2u+2*T ' 2-S 
1760  L3u*L3u'*S*T 
1270  HE XT  K 
1280  L2u*L2u  2 
1290  L3u«L3u!4 
1380  L 1  v*L2v*l.  3u»0 
1310  V I » V  c I  ) 
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Table  D-2.  (Coat'd)  Program  for  Equlspaced  Planar  Array 


1320  FOR  K« 1  TO  Hy 

1330  S*Wy<IO 

1340  T«$*C0S<.<K-.3>*Vi  > 

1330  S*$A2 

1360  Llv*Llv+T 

1370  L2u*L2u+2#T''2-$ 

1380  L3u«L3u*S*T 

1398  NEXT  K 

1400  L2u*L2v'2 

1410  L3v*L3u/4 

1420  Lluu*Llu«Llu 

1430  L2uv*L2u*L2u 

144C  Id»*l  powtr  <  1  >“Lluu''2 

1430  Au«r*g«pou*r  <  I  •■Cl£*Idt«1  pow«  r  •:  I  >  % 

14  60  V«r_pow4r  *Com  +  F3*L luu*L3u*L3u  +  F2t>*L2u”' F  1  <  :  :  r-*l  £y»  "♦LluV' 

1470  S i gm*pow*r < I > "SQR"  V*r  pou*r> 

1480  NEXT  I 

1490  PLOTTER  1$  “GRAPHICS" 

1300  GRAPHICS 

1510  SCALE  0,Num,-?,0 

1520  GRID  Num'10, 1 

1530  PENUP 

1340  FOR  I«0  TO  Hum 

1550  A*Id*»lpou*r< I > 

1360  IF  A> IE-7  THEN  1590 
1570  PENUP 
1380  GOTO  1600 
1590  PLOT  I.LGTCA) 

1600  NEXT  I 

1610  PENUP 

1620  FOP  1*0  TO  Hum 

1630  PLOT  I , LGT<Aw*r aj*pow*r ( I >  ' 

1640  NEXT  I 

1650  PENUP 

1660  FOR  1*1  TO  Hum 

1670  PLOT  I  ,  LGT  <  A*<*r  *g*po'J«r- <  I  >  ♦£  i  gm*pcwt'~ 1  I  ' 

1630  NEXT  I 

1690  PENUP 

1700  FOR  1*0  TO  Hum 

1710  PLOT  I  ,  LGT'Av*r  &g«pou€T  <  I  >t2*S  i  gm»p.iJ«ir  1  • 

1720  NEXT  I 
1 730  PENUP 
1740  PAUSE 
1730  DUMP  GRAPHICS 
1760  PRINT  L I N <. 8  > 

1770  PRINTER  IS  16 
1780  END 
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